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Abstract 

Using a recent technique, proposed by Eardley and Giddings, we extend 
their results to the high-energy colhsion of two beams of massless particles, 
i.e. of two finite-front shock waves. Closed (marginally) trapped surfaces 
can be determined analytically in several cases even for collisions at non- 
vanishing impact parameter in D > 4 space-time dimensions. We are able to 
confirm and extend earlier conjectures by Yurtsever, and to deal with arbitrary 
axisymmetric profiles, including an amusing case of "fractal" beams. We 
finally discuss some implications of our results in high-energy experiments 
and in cosmology. 
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I. INTRODUCTION 



The gravitational collapse induced, in classical General Relativity (GR), by trans- 
Planckian-energy collisions of particles and/or waves has attracted much theoretical at- 
tention since the early seventies. In the pioneering papers by Khan and Penrose [|l| and 
by Szekeres p| (see also 0, [Q), the case of infinitely-extended homogeneous plane waves 
was solved analytically in the interaction region where a (naked) singularity is inevitably 
produced. At the other extreme, the scattering of point-like objects (or of black holes) at 
zero impact parameter was investigated in classic papers by D'Eath and Payne [^, while 
R. Penrose managed to obtain a rigorous lower bound on the fraction of incident energy 
ending up in the black hole inevitably resulting from the collision. The difficult, intermediate 
case of the collision of finite-front shock waves has received comparatively little attention, a 
noteworthy exception being ref. 

Somewhat more recently, trans-Planckian scattering of particles and strings were investi- 
gated at the quantum level, as a gedanken experiment aimed at answering some fundamental 
questions in quantum gravity and/or in string theory. The problem turned out to be tractable 
either at large impact parameters (6 ^ -R^ ~ {GEy^^^~^^) through an eikonal approxi- 
mation, or when string-size effects manage to screen the non-linear classical effects that 
should trigger a collapse. In either case no black hole is formed. In spite of much effort, the 
region b < Rg, where black holes are expected to form, has remained untractable. 

The renewed recent interest in the field stems mainly from the following motivations: 



String-inspired cosmological models, such as the pre-big bang scenario |T^, connect 



(dilaton-driven) infiation to gravitational collapse through a conformal change of the 
metric. In order to avoid fine-tuning the initial conditions, it is crucial that collapse 
occur as generically as possible. The case of spherical symmetry was addressed in [p!T[] , 
while that of exact planar symmetry was solved analytically in 0], using precisely 
the techniques of 0, and 0. Since exact spherical or planar symmetry are quite 
special, it looks very desirable to extend the calculations to the case of finite-front 
shock waves, but, so far, very little progress was made, if any. 



The idea of large extra dimensions and of the brane Universe [|13[ allows for gravity 
to become higher-dimensional, and stronger than usual, below almost macroscopic 
distances. The true scale of quantum gravity could become as low as a few TeV. In 
such a context, black-hole formation in one of the near- future accelerators is all but 



excluded (for recent work on the subject, see and |T^ and references therein), but 



there has been some debate lit] as to the actual value of the cross section for black- hole 
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formation. 



New cosmological models based on the brane Universe idea have recently been pro- 
posed [l7]. In these models the big bang event is identified with the instant at which 



two almost parallel branes collide. Although the branes move slowly in this case, it 
is possible that techniques used for the relativistic case can be generalized to other 
situations endowed with similar symmetries. 

Recently, Eardley and Giddings [|18| proposed a promising technique for determining the 
occurrence of closed trapped surfaces (CTSs) in the collision of two shock waves. They 
considered the case of two colliding point particles at generic impact parameter b and in 
any D > 4. They succeeded in constructing CTSs for general b in D = 4, and at 6 = 
for D > 4, and offered arguments in favour of the existence of a CTS also in D > 4 for 
sufficiently small b. A consequence of their results is a lower limit on the cross section for 
BH formation from point-like particles. However, it is well known that strings behave 
rather as quasi-homogeneous beams over a size of order A^, the (quantum) string-length 
parameter. As we shall discuss at the end, this could reduce somewhat the lower bound on 
the cross section. 



With all these motivations in mind we shall extend the work of ||T8[ to the case of finite- 
size beam-beam collisions. We will first review, for completeness, the argument of [|l8l and 
then apply it to the case of homogeneous beams of finite size, first for D = 4 and 6 = 0, 
then for D = 4 and 6 7^ and then for D > 4 and any b. Finally, we shall present results for 
axisymmetric collisions when the profile of the beam is arbitrary and discuss some physical 
implications of our results. 



II. A CRITERION FOR GRAVITATIONAL COLLAPSE 

An impulsive wave moving with the speed of light along the positive z axis of a 
dimensional space-time leads to the metric (see e.g. 0, [pO[| ) 



ds = —dudv + (f){x.)6{u) du + rfx , (2-1) 

where 11 = 1 — z, v = t + z and x are the d = {D — 2) transverse coordinates. Einstein's 
equations require 

A(/)(x) = -167rGp(x), (2.2) 
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where p(x) is the energy density (energy per unit transverse hypersurface) in the beam. Let 
us use from now on units in which SttG = 1 and p has (for any D) dimensions of an inverse 
length. 

In the coordinates m, x the geodesies and their tangent vectors are discontinuous across 
the shock [Q, |20|, whereas in the new coordinates 



u = u 



1 



v = v + (t){^)e{u) + ^u^(m)(V0(x))2 



u 



x' = x' + -^(u)Vi0(x) 



(2.3) 
(2.4) 

(2.5) 



they are continuous. In these coordinates the metric becomes |]T8 

ds^ = —du dv + HikHjk dx^ dx^ , 



(2.6) 



where 



Hij = Sij + i\/ j(l){yi)u9{u) 



(2.7) 



Let us now consider the collision of two particle beams, or shock waves, moving in opposite 
directions along the z axis. By causality, outside the interaction region m > 0, f > 0, the 
metric is given by a trivial superposition of two metrics of the form 



ds = —du dv + 



dx^ dx^ , 



(2. 



where 



Hil'^ = + 7TViV,-0i(x)M^(«) 



l,2(x) = -2pi_2(x). 



(2.9) 

(2.10) 
(2.11) 



We reproduce now, for completeness, the construction of |T8[ to find a (marginally) closed 
trapped surface (CTS) S lying in the union of the two null hypersurfaces u = 0,v < and 
V = 0,u < 0. This property of S allows us to use the simple block- diagonal metric ( [2.81 ). 
On the null hypersurface u = 0,v < 0, the non-vanishing components of the Christoffel 
connection are simply 



F'^,^. = 2r\^ = V.V,0i(x). 



f2.12) 
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Actually, this result is strictly valid for n > e > with e arbitrarily small, so that we 
resolve a possible ambiguity by defining 6{0) = 6{e) = 1. Similarly, on the null hypersurface 
V = 0,u < 0, the non-vanishing components of the Christoffel connection are 



ViV,02(x) 



(2.13) 



Let us define S = Si U S2, with Si. u = 0,v = —■ipi{'x) < 0, and S2: v = 0,u = 
—ip2{'^) < 0; then S intersects the d-dimensional hypersurface m = = f on a closed {d — 1)- 
dimensional hypersurface C. We recall that a CTS is a closed space-like d-dimensional 
(hyper) surface S such that the two families of null geodesies orthogonal to S are converging 
at S. Modulo other conditions on the energy-momentum tensor that are met in our case, 
the existence of a CTS guarantees the occurrence of gravitational collapse, i.e., typically, the 
emergence of singularities hidden behind black-hole horizons. Rather than for CTSs we will 
look for marginally trapped (hyper) surfaces (MCTSs), on which the above null geodesies 
have zero convergence. In general S will be defined by /i = 0, /2 = 0, where fi and /2 are 
C^-functions such that /i;^ and f2;u are non-vanishing, non-parallel, and satisfy 







(2.14) 



for two distinct real values fii and fi2 of fi. Let A^f and ^^2* be two null vectors normal to S 
and proportional to g'"'{fi-u + f^if2;u) and g'^'ifi-u + I^2f2;u), normalized by N^N^g^^, = -1 
and let F^f (a = 1,2, ... ,d) be a set of space-like unit vectors orthogonal to each other and 
to and N2. The two null second fundamental forms of S are defined by 



Xn^u 



(2.15) 



iS is a CTS (MCTS) if g^'^Xifj.u and g'^^X2^iu are never positive (vanish) on S. On Si we 
choose 



/ \ 

-1 




V I 



( 



-2 \ 



2 

^1,1 



( \ 

-^1,1 
1 


V I 



5 2 



One then easily finds that 

g'"'xi,u = Q, /'^X2M. = A(0i-v^i). 



/ \ 

-V'1,2 

1 

I ; 



(2.16) 



(2.17) 
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On ^2 we choose, analogously. 







\ I 



2 ' 

-2 

^2,2 



\ ^2,d J 





1 



V / 



( -V'2,2 \ 




1 

V I 



giving 



Continuity of the outer null normal N2 on C (u = f = 0) requires 

_ V^2 _ 
V'2,1 ^2,2 

i.e. Vipi to be parallel to V'?/'2 and 

{V^if{V^2f = 16. 



The necessary and sufficient condition for this to happen is 

V^i ■ Vi)2 = 4. 

To summarize, iS is a MOTS under the following conditions: 

ipi^2 > inside C, ipi^2 = on C 

^(V'1,2 — 01,2) = inside C 
Vipi ■ V^2 = 4 on C. 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



(2.23) 
(2.24) 
(2.25) 



In [IS] the authors construct C for the collision of massless point particles at zero impact 
parameter b and any D, as well as for non- vanishing impact parameter in D = 4. In the 
following sections we will describe extensions to the case of non-point-like beams. 



III. HOMOGENEOUS FINITE-SIZE BEAMS 



Let us consider the case of finite-size beams with radius Ri > R2 and homogeneous 
energy density pi 2 inside. It is useful in this case to introduce for each beam its "focal 
distance" 
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1,2 



2^1,2, 



(3.1) 



This is where the null geodesies parallel to the z-axis converge after hitting the shock wave, 



which therefore acts as a perfect anastigmatic lens [21], |20| 



A. D = 4, 6 = 

For 6 = 0, = 4 ((i = 2) we choose 



2/i,2V'i,2(r) = (i??,2 - r')^(i?i,2 - r) - log (^-^ j e{r - i?i,2) + Ci,2 (3.2) 
AV^i,2(r) = -2pi,2^(i?i,2-r), (3.3) 



where r is the radial coordinate in the transverse x space. The conditions (|2.24| ) are already 
satisfied while (|2.25|) reads 



4 = ^Mr)^Mr) = -^0{R2 - r) + -^e{R, - r)e{r - R,) + ^^e{r - R,). (3.4) 
If 



i?2>2v//i/2, (3.5) 

( p.25| ) has two solutions: 

r = r,i = 2/A^<i?2 , (3.6) 

r = r,2 = > Ri. (3.7) 

The first lies inside both beams, while the second is external to both. The physical meaning 



of this result is quite clear: if R2 > 2^^/1/2 there are CTSs that intersect the collision plane 
at any value of r between r^i and (the CTS becomes a MCTS at r = r^, rc2). An example 
is the surface S = SiU S2, with 

S^:u = v + {rl- r^)lr^ = , 2fi < r < Ri, (3.8) 

so that g'^''X2nv = 4(r^^ — (2/i)^^) < and similarly for ^2. Thus 5 is a CTS for 
2max(/i, /2) < rc < R2. 

In the case 
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R2 — /1/2 , 



(3.9) 



we have MCTSs intersecting the colhsion plane for any r in i?i > r > i?2- The conditions 
( p.23| ) finally determine the constants Ci^2- So we have shown that in the collision of finite-size 
beams with vanishing impact parameter in 4 dimensions a CTS forms if R2 > 2^^/1/2. This 
result is in full agreement with one of the conjectures by Yurtsever [0. In fact his criterion 
for collapse, R2 ^ 2a//i/2, is simply replaced by R2 > c^/flf2 with c < 2. 



B. L> = 4, 6/0 

For the case of non-vanishing impact parameter, a closed trapped surface can be con- 
stucted starting with the solution at 6 = 0. Let the first beam be centred around —a < 
on the X-axis and the second beam around a on the x-axis, so that the impact parameter is 
b = 2a. 

Let us construct a solution starting with ( p.6|) . If the impact parameter fulfils the condi- 
tion 



&<2i?2- 4^/1/2, (3.10) 

then on the circle 

C:x2 + y2 = 4/1/2 (3.11) 

we are still inside both beams. This time (unlike in the case 6 = 0) we can allow in ^/^i 2 — 0i,2 
not only a constant but also a term linear in x. We thus choose 

2/i,2^i,2(x, y) = -{{x ± a)2 + - R^eiRl^ - {x ± - y^) 

-Rh log ^^^#^^((^ ± + y'- RI2) 
^1,2 

+ [4/1/2 - RI2 ± 2aa; + a^] , (3.12) 

so that 

AV^i,2(x, y) = -2pi,2^^(i??,2 -ix± af - y^) (3.13) 

V'i,2(x,y) = 0onC (3.14) 

Wi ■ V7/'2 = 4 on C. (3.15) 



8 



We conclude that if (|3.5| ) and ( p.lOj ) hold, a black hole will form in the collision, the maximal 
impact parameter being 



bmaxA = 2-R2 — 4a//i/2 



(3.16) 



A similar method can be used to generalize ( |3.7] ). Let us assume, for simplicity, that the 
two beams are identical: Ri = R2 = R, Pi = P2 = P- If the impact parameter fulfils 



p2 

< 6 < y - 2i?, 



(3.17) 



then the circle 



C -.x^ + 



R^ 
4/2 



(3.18) 



is external to both beams. Choosing 



2M,2(a;, y) = -{{x ± + y^- R^iR^ - {x ± af - y^) 



-i?^logi"^")' + ^' 



+R^ log 



e{{x ± af + y^-R^ 

2„,2 Jl 



% IR^ 



• c2 



(3.19) 



we find 



AV^i,2(a;,y) 
^i,2(x,y) 



-2pQ{R^ 
on C 



(x ± af - y') + pnR'dix ± — j5(y) 



V^i ■ V^2 = 4 



(r?2 - a^) 



^c2 



on C. 



(3.20) 
(3.21) 

(3.22) 



The extra sources on the x-axes at always lie outside C. The situation is now the same 
as for the collision of two particles described in [|1^ . We use complex variables , 



x+iy 



and 



make a conformal transformation 



z\z)= log - 



2— 



1 - 

rc2 



(3.23) 



The circle C is then transformed to a new curve C and, in terms of the new coordinates. 



9 



AVi,2(a;',?/') = -2pB{R^ - {x' ± af - y'^) inside C (3.24) 

^^ 2(a;',i/') = on C (3.25) 

VVi ■ Vi)2 = 4 on C. (3.26) 

The points (±a, 0) are transformed to (±a', 0) = (±|rc2/i(a/rc2), 0) where 

1 — w"^ 1 + w"^ 

h{w) = ^log-^ (3.27) 

W 1 — 

has a maximum at Wmax ~ 0.62519 with h{wmax) ~ 0.80474. Taking (|3.17|) into account we 
have that this solution is possible for 

U -0 ' <^ U -I ^^^2 h{w^ax) for R<{1- Wmax) ^^2 

- ~ I re2 Ml - R/rc2) for i? > (1 - w^^^) r,2 " ^"'■'''^ 
Note that (|0| ) only fixes i? < rc2 and that 6rraax,2 > bmax,i- 



C. D > i, b = 

Let us generalize the results to arbitrary dimensions D > 4. For vanishing impact 
parameter we find (recalling the definition / = {d/2){87iGp)^^): 

2/i,2^i,2(r) = -(r^ - Rl,)eiR,,2 - r) + 2R^^^^ ''' 0{r - Ri,2) + Ci,2 (3.29) 

AV'i,2(r) = -2pi,2e(i?i,2 - r). (3.30) 



The condition (|2.25|) reads 



4/1/2 = -^Mr)-^Mr)fif2 = r'e{R2 - r) + r'-^R^-^e{R, - r)e{r - R^) 

+ r6-2^i?f-2i?f "2^(r - i?i). (3.31) 

Keeping in mind that, by convention, R2 < Ri, let us introduce two length scales: 

f = 4hh,F = fiy-£^ >/. (3.32) 

We then distinguish various cases according to the value of i?2- 

• R2<2f 

In this case we cannot find any MCTSs. 
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2/ < i?2 < 2F 

We find two MCTSs: the first intersects the colhsion hyperplane at 

r = r,i = 2/ < i?2, (3.33) 

i.e. inside the two beams. The second intersects at 

- = = j ' (3-34) 

i.e. inside beam 1 and outside beam 2. 
i?2 > 2F 

We find again two MCTSs: the one corresponding to the circle of radius Vd and a 
second with 

r = r.3=(^^i0^)'^>i?„ (3.35) 
i.e. outside both beams. 



Finally, the constants ci,2 can be easily determined from ( |2.23| ). 

We have thus shown that in the collision of finite-size beams at vanishing impact pa- 
rameter in {D > 4) dimensions, a closed trapped surface will form whenever i?2 > 2/. The 
radii rd and rc2,3 all correspond to MCTS and, by continuity, there should be genuine CTS 
intersecting the collision hyperplane at r^i < r < rc2,3. 

It is amusing to try to understand, at any D, the physical meaning of the critical radius 
which lies inside both beams. Boost the system to a Lorentz frame in which pi = P2 = P 
and go back for a moment to the coordinates u, v and x to analyse the null geodesies. A 
null geodesies parallel to the 2;-axis and initially at m < 0, f = and R > xq > will hit the 
shock wave at m = 0, and then instantaneously jump by 



At = Az = ^ = (3.36) 

and again follow a straight line to reach a; = at 

vp = 0, tF = -zp = ^ = f- (3.37) 

Precisely for xq = 2/ we have Az = zp and therefore a scattering angle of 7r/2 in the 
(a;-2;)-plane. In other words, for R2 = 2/ the lens effect of the wave is strong enough to 
deflect the energy impinging on its edges by 90 degrees! Fig.0 illustrates this point. 
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T>. D > i, b^O 



For non-vanishing impact parameter, the solution starting with (|3.33| ) is constructed as 
in 4 dimensions. If the impact parameter fulfils 



(3.38) 



then the circle 



C:xl + --- + xl = 4/1/2 



(3.39) 



lies inside both beams and we can take 



- -((xi ±af + xl + ■ ■■ + xl- RI2) 
x\ — ■■ ■ — x^) + 9{{xi ±af + x\ + h X, 



-^1,2) 



[{xi ± a)^ + X2 + h x^ 



4-D 
2 



A-D 



+ 



4/i/2-i??2±2axi + a2 



(3.40) 



in order to satisfy all the conditions for a MCTS. 



Therefore if R2 > 2/ and ( p.38| ) hold, a black hole will be created in the collision. To 
find a solution starting with (|3.34|) is more difficult, but it is easy to argue, by continuity, 
that an outer MCTS should exist also in this case. 



IV. NON-HOMOGENEOUS BEAMS 



We now want to address the question of how natural the creation of a black hole is in the 
collision of two particle beams of arbitrary transverse profile. We will restrict our attention 
to axisymmetric examples, i.e. to collisions at 6 = of beams whose energy density is only 
a function of the transverse radius pi^2(x) = pi,2(^)- Then the beam energy inside a radius 
r is simply 

r 

E{r) = Qd~2 Jdr' r'^-^p{r'), (4.1) 


where ild = 27r'^/^/r((i/2) is the solid angle in d dimensions. We want to construct a closed 
trapped surface and for this we start with 

Mr)=Mr) + Cu Mr) = Mr) + C2, (4.2) 
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so that (|2.24|) is fulfilled. This also implies 



so that ( |2.25|) can be written as 



1 = o \D-3 ^E,ir,)E2ir.)- (4-4) 

If this has at least one solution, we can adjust the constants ci^2 in such a way that ?/'i,2 
vanishes on the critical radius, i.e. we have constructed a MOTS. 

As a first example, consider a "fractal" beam 

Pi,2(r) = ai,2r'^-^(5>0, (4.5) 

where ai_2 are constants and the last condition guarantees that E{r) ~ < oo for finite r. 

We can call 6 the fractal dimension of the (energy stored in the) beam. For 5 7^ — 1 we 
always find a critical radius for a MCTS at 

52 



S-d+l 



(4.6) 



while, for 5 = d — 1, we can choose any radius if ai^2 fulfil the condition 

v/^=(rf-l)^ (4.7) 

but we find no solution otherwise. We interpret this by saying that, for 6 > d — 1 {6 < d—1) 
we have CTSs with r > Tc (r < Tc), with Tc given by (|4.6| ), while, for 6 = d — 1, we either 
have, at all r, a MCTS (if (g^) holds), a CTS (if > {d - if) or, finally, nothing at 

all, if ^/ala2 < {d — 1)^. 

Consider finally the case 

p(r) = ai,2 r^-V^i'^'^', (4.8) 

where ( |4.4| ) becomes 



10:10:2 



1 -e-^i '■')(1 '^?), (4.9) 



and the existence of a solution depends on the precise values of 01^2, ^1,2 and D. As an 
example, for D = 5, 61 = 62 = S, oi = 02 = a = 26 there is a solution Tc > for 5 > 1 but 
none otherwise. 
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V. PHYSICAL IMPLICATIONS 



Let us briefly discuss the possible physical implications of our results, first in string 
cosmology and then on cross-sections for black-hole formation. 

• String cosmology 



We believe that our results clearly go in favour of the arguments given in |]TT| , |]T2[ for 
the generic occurrence of gravitational collapse in the presence of an initial, classical 
chaotic sea of massless waves. It is also self-evident that the criteria for collapse never 
involve the Planck length/mass, but only dimensionless ratios of classical lengths de- 
scribing either the geometry of the beams or the geometry of space-time. We have thus 
provided further answers to the allegations of fine-tuning in pre-big bang cosmology 
claimed in One limitation of our method is the restriction to impulsive waves. 

We think that this should not be a problem of principle and that generalization of our 
results to "thick" waves should be possible. By contrast, our results have nothing to 
say on the nature of the singularity that lies inside the CTSs and, in particular, on 
whether, in its vicinity, space-time is described by a Kasner-like metric or by the more 
generic BKL oscillatory behaviour recently discussed by Damour and Henneaux [^. 
This may very well depend on the nature of the collapsing waves. 

Black-hole production at accelerators 



Our results confirm the absence of the exponential suppression claimed in [|T6[ in agree- 



ment with the original estimates in |14] and more recent work |25] and \2€.]. However 
they imply a revision, unfortunately downwards, of previous estimates and [|r^ of 
black hole production in theories with large extra dimensions and low-scale quantum- 
gravity. Those theories make sense only in a superstring context, which introduces a 
length scale A^, which is at least as large as the (true, multidimensional) Planck length 
Ip ~ Mp^ . As mentioned in the introduction, if strings, rather than point particles, are 
colliding, string-size effects can be modelled by considering beam-beam collisions 
with beam sizes of order As. Using our results ( |3.1CI| ), ( |3.38| ), and inserting i? ~ As, we 
find that the range of impact parameter where black hole formation must occur is 

|6| < 0(i?s - As) , Rs ~ (GE)V(^-3), (5^1) 

where E is the centre-of-mass energy of the collision. Equation ( |5.1D , together with the 
expected relation |T5| ubh ~ ^rfo^, implies that, in order to arrive at reasonable cross 



sections for black hole production, the cm. energy should satisfy E > -^threshold ~ 
Mp{Xs/lp)^^^, i.e. should be parametrically larger than Mp if, as expected, As//p > 1. 
This is in agreement with many of the conclusions reached in (see also [p7[] ). 
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FIG. 1. Null geodesies in the metric ( |2.lD for homogeneous beams with /i = /2 = / = -f (from 
Ref. [^) and the hypersurface 5i C 5. The geodesies come in from the left, parallel to the z-axes, 
then jump according to the dotted lines and reappear on the circle (j + 1)'^ + (j)^ = 1. All the 
geodesies hitting the wave at t = z = converge in — / at the same time, the outermost coming in 

2 

at 90 degrees. The bold line is the hypersurface Si : z = t = jj — f . The dashed-dotted line is 

2 

the hypersurface z = —fr. 
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